arXiv:1508.01728vl [astro-ph.HE] 7 Aug 2015 


Possible identifications of newly observed magnetar 
quasi-periodic oscillations as crustal shear inodes 


Hajime Sotani 1 , Kei Iida 2 , Kazuhiro Oyamatsu 3 

1 Division of Theoretical Astronomy, National Astronomical Observatory of Japan, 2-21-1 

Osawa, Mitaka, Tokyo 181-8588, Japan 

2 Department of Natural Science, Kochi University, 2-5-1 Akebono-cho, Kochi 780-8520, 

Japan 

3 Department of Human Informatics, Aichi Shukutoku University, 2-9 Katahira, 
Nagakute, Aichi 480-1197, Japan 


Abstract 

Quasi-periodic oscillations (QPOs) discovered in soft-gamma repeaters (SGRs) 
are expected to help us to study the properties of matter in neutron stars. In 
earlier investigations, we identified the QPOs of frequencies below ~ 100 Hz 
observed in giant flares of SGR 1806—20 and SGR 1900+14 as the crustal 
torsional oscillations. For this purpose, we calculated the frequencies of the 
fundamental torsional oscillations with various angular indices £, by changing 
the stellar mass and radius. In this work, we try to explain the additional 
QPO frequencies recently reported by Huppenkothen et al. 0,3 within the 
same framework as before except that we newly take into account the effect of 
electron screening, which acts to decrease the frequencies by a small amount. 
Those QPOs were discovered in two different SGRs, i.e., SGR 1806—20 and 
SGR J1550—5418. Then, we find that the newly observed QPO frequency 
in SGR 1806—20 can be still identified as one of the frequencies of the fun¬ 
damental torsional oscillations, while those in SGR J1550—5418 can also be 
explained in terms of the torsional oscillations although the relevant angular 
indices are difficult to identify. 
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1. Introduction 


Neutron stars, which are stellar remnants of the core collapse of massive 
stars, give us one of the best opportunities to investigate the physics under 
extreme conditions. This is partly because the density of matter inside the 
star can become significantly larger than normal nuclear density, partly be¬ 
cause some neutron stars are strongly magnetized while others are rapidly 
rotating, and partly because the gravitational fields around neutron stars are 
strong enough for effects of general relativity to manifest themselves. Astro¬ 
nomical phenomena associated with neutron stars could leave imprint of the 
physics under such extreme conditions. The asteroseismology is a powerful 
technique to reveal the interiors of neutron stars via their oscillation spectra, 
as in the case of the seismology in the Earth and the hclioseismology in the 
Sun. In fact, through detection of the frequencies of neutron star oscillations, 
it might be possible to deduce neutron star masses and radii, the equation 
of state (EOS) of matter in the star, and so on (e.g., BBS 00). Gravi¬ 
tational waves, if detected, could be one of the most promising sources that 
provide the oscillation spectra of neutron stars. 

Alternatively, it is expected that observational evidences for neutron star 
oscillations have already been given as the quasi periodic oscillations (QPOs) 
observed from soft gamma repeaters (SGRs) in the afterglow of the giant 
flares. Up to now, at least three giant flares were observed, and various QPO 
frequencies were found in two of them, which radiated from SGR 1806—20 
and SGR 1900+14 0 0 IQ]- Since SGRs are considered to be magnetars, 
which are strongly magnetized neutron stars, the observed QPOs are ex¬ 
pected to be strongly associated with the neutron star oscillations. There 
are many attempts to explain these QPOs theoretically, which are based on 
shear torsional oscillations in the crustal region of a_ neutron star and/or 
magnetic oscillations throughout the star (e.g., [ll|, 12, EiHEBEQ)- 
Subsequently, analyses of elastic-magnetic oscillations in magnetars reveal 
that, depending on the magnetic field strength, the oscillations near the 
star’s surface are excited basically by either crustal torsional oscillations or 
magnetic oscillations [l8|, 13- According to observational estimates of the 
surface magnetic field strength of the SGRs from which the giant flares radi¬ 
ate [ 20 I, 0 , the magnetic field in the star can give rise to a restoring force, 
whose strength is comparable to that induced by the crustal shear modulus, 
and hence shear torsional oscillations and magnetic oscillations are almost 
indistinguishable. However, given that the highest magnetic fields are more 
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or less localized, e.g., in a toroidal form, it is reasonable to start with the 
assumption that the QPOs observed in SGRs are identified as the crustal 
torsional oscillations. Then, such an identification could tell _us_ information 
about the properties of the crust, particularly the EOS [ 22 ], 


23, 24 


One of the most important parameters characterizing the EOS of neutron- 
rich matter in the crust is the density dependence of the nuclear symmetry 
energy, L , which is strongly associated with the thickness of the region where 
nuclei exist as non-uniform nuclear structures [ 25 ], 26, 27 


The constraint 
1, but it is still 


on L can be given via the terrestrial nuclear experiments [2 
difficult to obtain a severe constraint. On the other hand, we developed a way 
of constraining L through the identification of the QPOs observed in SGRs 


as the crustal torsional oscillations [29|, |30|, [31] as well as through possible 


simultaneous mass and radius measurements of low-mass neutron stars [32]. 
This is a constraint from nuclear matter at extremely large neutron excess, in 
contrast with the case of the terrestrial nuclear experiments associated with 
nuclear matter at relatively small neutron excess. 

Recently, new QPO frequencies in SGR J1550—5418 and SGR 1806—20 
have been reported. In SGR J1550—5418, the QPO frequencies of 93 and 
127 Hz have been discovered from a storm of 286 bursts [1]. In addition, the 
QPO frequency of 57 Hz has been discovered in the shorter and less energetic 
recurrent 30 bursts radiating from SGR 1806—20 [ 2 ]. Since the QPO of 
frequency 57 Hz comes from the same SGR as that from which the giant 
flare was detected, it is important to make sure that this new QPO frequency 
can be explained within the framework that reproduces the low-lying QPO 
frequencies found in the giant flare. Meanwhile, although SGR J1550—5418 
adds to a list of the SGRs that have the QPOs detected, unfortunately the 
observed QPO frequencies are limited and not low enough to clearly identify 
the corresponding angular indices. Thus, we here systematically examine 
how valid the crustal torsional oscillations are to explain the QPO frequencies 
observed in the SGRs. In order to calculate the frequencies of the crustal 
torsional oscillations, we adopt the shear modulus that allows for the effect 
of electron screening on nuclei in the crust [33j. Additionally, as the effect of 
neutron superfluidity on the enthalpy density of matter in the inner crust, we 
adopt the results for the superfluid density given by Charnel [34jJ, as in Refs. 
30l . [31]. We adopt the geometric unit of c = G — 1 in this paper, where c 


and G denote the speed of light and the gravitational constant, respectively. 


3 

























2. Crust configuration and EOS parameters 

In the vicinity of the saturation point of symmetric nuclear matter at 
zero temperature, the bulk energy of nuclear matter per nucleon w can be 
generally expressed as a function of baryon number density n b and neutron 
excess a, as in Ref. [35]: 


w = w 0 + 


Kp 

1871q 


(n h - n 0 ) 2 + 



L 

3 no 


(n h - no) 


a 


( 1 ) 


where w 0 , n 0 , and K 0 are the saturation energy, saturation density, and 
incompressibility of symmetric nuclear matter (a = 0). In addition, So 
and L are associated with the density dependent symmetry energy S(n^) 
as So = S(no) and L = 3n 0 (dS/dn h ) nh=no . Note that only wo, no, and 
So among the five parameters are well constrained from empirical data for 
masses and radii of stable nuclei [36]. 

In describing matter in the crust of a neutron star, we follow the derivation 
by two of the authors (K.O. and K.I.), which is based on a phenomenological 
approach j27j[. Hereafter, we refer to the resultant phenomenological EOS 
of matter in the crust as OI-EOS. The OI-EOS is constructed as follows. 
First, various models for the bulk energy w{n-o,a ) of nuclear matter are 
made in such a way as to reproduce Eq. (JT]) in the limit of no —> 0 and 
a —> 0. Then, within a simplified version of the extended Thomas-Fcrmi 
theory, the density profile of stable nuclei was obtained for each model for 
w. Finally, the optimal values of u>o, no, and So were determined so that the 
charge number, mass excess, and charge radius that can be calculated from 
the density profile obtained for given y = —KoSo/{3n 0 L) and K 0 should 
fit well to the experimental data [36]. After that, in order to obtain the 
equilibrium nuclear shape and size as well as the crust EOS for various sets 
of y and Kq, we generalized the Thomas-Fcrmi model by adding dripped 



For given mass M and radius R of a nonrotating neutron star, the crust 
configuration can be constructed by integrating the Tolman-Oppenheimer- 
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Volkoff equations with the above-mentioned crust EOS from the star’s sur¬ 


face inward up to the basis of the crust [37], [29j, l3(J, l31| . This construction 


effectively avoids uncertainties in the core EOS. In this work, we particu¬ 
larly consider typical neutron star models with M and R in the range of 
1.4 < M/M© < 1.8 and 10 km < R < 14 km. 

One of the crucial properties that govern the shear oscillations is the 
elasticity, which is characterized by the shear modulus n . The shear modulus 
in the crust is mainly determined by the lattice energy due to the Coulomb 
interaction, which is approximately given as 


H = 0.1194 x 


rii(Ze) 


( 2 ) 


where n*, Z, and a = (3/47Trq) 1 / 3 are the ion number density, the nuclear 


charge number, and the radius of a Wigner-Seitz cell, respectively 381 ]. This 


formula is derived in the limit of zero temperature from Monte Carlo cal¬ 
culations averaged over all directions on the assumption that the nuclei are 
point particles forming the body-center cubic lattice [39]. 

On the other hand, Kobyakov and Pethick proposed the modification of 
the shear modulus, where the effect of electron screening is taken into account 


33]. The modified shear modulus can be expressed as 


/i = 0.1194 (1 - 0.010Z 2/3 ) 


rii(Ze) 


( 3 ) 


where the term of Z 2 / 3 represents the effect of electron screening. That is, 
the shear modulus decreases due to such additional effect, which in turn leads 
to reduction of the frequencies of torsional oscillations 40 Since the charge 
number Z also depends on L, however, it is still unclear how the reduction of 
the frequencies of torsional oscillations due to the effect of electron screening 
depends on L. At subnuclear densities, Z tends to decrease with L, because 
the smaller symmetry energy, corresponding to larger L, helps more protons 
to change into neutrons. In fact, Z decreases with L for the OI-EOS adopted 
in the present work, as shown in Fig. [0 Consequently, one expects that the 
frequencies of torsional oscillations by using the modified shear modulus [Eq. 


1 The electron screening acts to change the toroidal oscillation frequencies also via mod¬ 

ifications of the enthalpy density and the crustal structure, but such change is negligibly 

small. 
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Figure 1: (Color online) Nuclear charge number, Z, plotted as a function of the energy 
density, p, for the OI-EOS with various values of L. The labels on the solid lines denote 
the values of L; L = 52.2, 73.4, 97.5, and 146.1 MeV correspond to the OI-EOS with 
(y,K o )={-220 MeV fm 3 , 180 MeV), (-220 MeV fm 3 , 230 MeV), (-220 MeV fm 3 , 280 
MeV), and (—220 MeV fm 3 , 360 MeV). For reference, the result from the EOS proposed 
by Douchin & Haensel 41 is also shown with the broken line, where the corresponding 
value of L is 45.9 MeV. 


0] differs little from those by using the original shear modulus [Eq. (121)]. 
if L is sufficiently large. Even so, we will examine the torsional oscillations 
with the modified shear modulus [Eq. ([3])] in this work. 


3. Crustal torsional oscillations in neutron stars 


We proceed to calculate the fundamental frequencies of the torsional os¬ 
cillations that are excited in the crust of a nonrotating neutron star of mass 
M and radius R as constructed from the crust EOS of a given set of L and 
K 0 in the previous section. In particular, we adopt the relativistic Cowling 
approximation to calculate the frequencies. That is, we neglect the metric 
perturbations and keep them zero during the oscillations. This is presumably 
a good approximation for considering the torsional oscillations, because they 
are axial parity oscillations, which do not involve the density perturbations. 
The perturbation equation for the torsional oscillations can be obtained from 


the linearized relativistic equation of motion 42]. Then, imposing the appro¬ 


priate boundary conditions at the top and bottom of the crust, the problem 
to solve reduces to the eigenvalue problem, where the eigenvalues correspond 
to the eigen-frequencies of the torsional oscillations. The perturbation equa¬ 
tion and the boundary conditions are explicitly described in Ref. jilj; pasta 


6 










nuclei, if present, are assumed to have zero shear modulus. 

Additionally, the effect of dripped neutrons should be taken into account 
for the calculation of the frequencies of the torsional oscillations. In fact, 
neutrons are generally considered to drip out of nuclei when the density 
becomes higher than ~4x 10 11 g/cm 3 , and some of them behave as a super¬ 
fluid. Unfortunately, it is still unclear how many dripped neutrons behave 
as such, but most of them can move together with the nuclei because of the 
entrainment effects 34]. In fact, according to the results from the band cal¬ 
culations in Ref. (34 1, only of the order of 10 — 30 % of the dripped neutrons 
can participate in superfluidity at n b ~ 0.01 — 0.4no- Through the enthalpy 
density, the frequencies of the torsional oscillations depends strongly on this 


ratio 30]; we adopt the results by Charnel 34] in the present work. Then, 
we calculate the frequencies based on the prescription how to build the effect 
of superfluidity into the perturbation equations shown in Ref. 


31 


First, to see the dependence of the frequencies on the parameters that 
characterize the crust EOS, L and K 0 , we calculate the frequencies of the 
fundamental torsional oscillations with l = 2 for the stellar models with 
M = 1.4M 0 and R — 12 km. The calculated frequencies are plotted as a 
function of L for various values of Kq in Fig. [2] Note that the shear torsional 
oscillations are often referred to as f-modes, which are labelled as n t(, with the 
angular index i and the number n of radial nodes in the eigenfunction. From 
this figure, one can observe that the frequency 0 t 2 is almost independent of 
the incompressibility K 0l while depending strongly on L. This tendency has 
been already shown in the case in which is calculated from the shear 

~30l ]. but still holds in the case of the modified 


modulus [Eq. ([2])] in Ref. 
shear modulus [Eq. 


Since we confirm that this is true for different 
stellar models within 1.4 < M/M 0 < 1.8 and 10 km < R < 14 km, 0 t 2 can 
be approximately expressed as 


0 t 2 — c 2 


(o) 


+ ef L\ 


(4) 

where c 2 °\ c£\ and cf' 1 denote positive coefficients that depend on the stel¬ 
lar models. In practice, this fitting formula agrees well with the calculated 
frequencies for various sets of the EOS parameters within the accuracy of a 
few per cent, as shown in Table [H 

In a similar way, we find that the frequencies oU for ^ > 2 are almost 
independent of Kq and can be approximately expressed as a function of L 
via 

(5) 


o ti = cf ^ — c» ±; L + cY’L 


(i) 


42) T 2 
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Figure 2: (Color online) The frequency of the fundamental torsional oscillations with i = 2, 
0 ^ 2 1 plotted as a function of L for the stellar models with M = IAMq and R = 12 km. 
The thick solid line denotes the fitting formula [Eq. ([1)1]. 


Table 1: The numerically calculated values of the frequency of the l = 2 fundamental 
torsional oscillations, o4°\ and the values obtained from Eq. a, o4 e) , for the stellar 
models of M = IAMq and R = 12 km. The relative errors determined by (ot^ — 
o4 e) )/o4 c) are also tabulated. 


y (MeV fm 3 ) 

K 0 (MeV) 

L (MeV) 

o4 c) (Hz) 

o4 e) (Hz) 

relative error (%) 

-220 

180 

52.2 

21.59 

21.12 

2.14 

-220 

230 

73.4 

17.79 

17.49 

1.71 

-220 

280 

97.5 

14.28 

14.21 

0.43 

-220 

360 

146.1 

10.46 

10.47 

-0.15 

-350 

180 

31.0 

25.25 

25.49 

-0.98 

-350 

230 

42.6 

22.89 

23.02 

-0.60 

-350 

280 

54.9 

20.45 

20.63 

-0.86 

-350 

360 

76.4 

16.65 

17.02 

-2.21 

-1800 

180 

5.7 

31.65 

31.68 

-0.12 

-1800 

230 

7.6 

31.30 

31.17 

0.41 

-1800 

360 

12.8 

29.87 

29.84 

0.08 
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Figure 3: (Color online) The frequencies of the fundamental torsional oscillations with 
l = 2, 3, and 4, plotted as a function of L for the stellar models with M = IAMq and 
R = 12 km. The broken lines correspond to the frequencies calculated from the shear 
modulus [Eq. lj2])]. while the solid lines correspond to those from Eq. lj3|). 

where cf\ c^\ and are adjustable parameters. 

In Fig.[3l we show the comparison between the frequencies calculated from 
the shear modulus [Eq. (j2J)] and those from the modified shear modulus [Eq. 
P)] for the stellar models with M = 1.4M 0 and R — 12 km. As discussed 
in the previous section, the deviations in o ti due to the effect of electron 
screening increase with decreasing L, but almost vanish for large L. 


4. Comparison with the QPO frequencies 


Let us now compare the calculated frequencies 0 tg with the QPO frequen¬ 
cies observed in SGRs. The identification of the QPO frequencies observed 
in SGR 1806—20 as the frequencies o tg is more difficult than that in SGR 
1900+14, because not only many QPO frequencies are discovered in SGR 
1806—20 but also the interval between the observed QPO frequencies 26 and 
30 Hz is remarkably small [Til ]. Nevertheless, as shown in Refs. jdol . 31 


the QPO frequencies 18, 26, 30, and 92.5 Hz observed in SGR 1806—20 
can be explained well in terms of the fundamental torsional oscillations with 
i = 3,4,5, and 15, where the frequencies of the torsional oscillations are 
calculated from the original shear modulus [Eq. (J2J)]. Here, we re-calculate 
the frequencies from the modified shear modulus [Eq. (|3|) ]. Then, we find 
that the same correspondence as in Refs. [3()|, 31] is still intact, as shown in 
Fig. [4] for the stellar models with M = 1.4M 0 and R = 12 km. We remark 
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Figure 4: (Color online) Correspondence between the QPO frequencies observed in SCR 
1806—20 and the calculated frequencies of the fundamental torsional oscillations for the 
stellar models with M = 1.4M 0 and R = 12 km. The horizontal broken lines denote the 
QPO frequencies originally reported by Israel et al. Strohmayer & Watts [Io|], while 
the horizontal dotted line denotes the new QPO frequency discovered by Huppenkothen et 
al. 0. The vertical line denotes the value of L = 123 MeV, with which the observed QPO 
frequencies agree best with the calculated frequencies of the crustal torsional oscillations. 

that the reproduction of the QPO frequencies by the crustal torsional oscil¬ 
lations is equally good. Moreover, we also include the new QPO frequency 
57 Hz reported by [2] in this figure, and we find that this additional QPO 
frequency can be identified as the £ = 9 fundamental torsional oscillations. 
As can be seen from Fig. [4j the most suitable value of L to adjust the crustal 
torsional oscillations to the five observed QPO frequencies is L = 123 MeV, 
with which the frequencies obtained from Eq. (J5]) and their relative errors 
from the observed QPO frequencies are shown in Table [21 Since not only 
the QPO frequencies discovered from the giant flare but also that observed 
from the different event in the same object, SGR 1806—20, can be explained 
in the same framework, we conclude that the torsional oscillations are still 
promising as the origin of the QPOs from SGRs. 

The same identification between the QPO frequencies observed in SGR 
1806—20 and the crustal torsional oscillations is also possible for different 
stellar models. In Fig. [5l we plot the optimal values of L to explain the 
observed QPO frequencies for various stellar models within 1.4 < M/Mq < 
1.8 and 10 km < R < 14 km. 

Unlike SGR 1806—20, the observed QPO frequencies in SGR 1900+14 
can be simply explained in terms of the crustal torsional oscillations. As 
in Ref. [30'], we End that the observed QPOs of frequencies 28, 54, and 84 
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Table 2: The QPO frequencies observed in SGR 1806—20, the corresponding angular 
indices and frequencies of the fundamental torsional oscillations, the latter of which are 
obtained from Eq. m by substituting the optimal value into L for M = 1.4Mq, and 
R = 12 km, and their relative errors from the observed QPO frequencies. 


QPO frequency (Hz) 

e 

o4 e) (Hz) relative error (%) 

18 

3 

18.62 

-3.44 

26 

4 

24.98 

3.92 

30 

5 

31.15 

-3.85 

57 

9 

55.22 

3.13 

92.5 

15 

90.76 

1.88 



Figure 5: (Color online) The optimal values of L to explain the QPO frequencies observed 
in SGR 1806—20, as in Fig. U for various stellar models. In the figure, the circles, dia¬ 
monds, and squares correspond to the results for the stellar models with R = 10, 12, and 
14 km. 
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Figure 6: (Color online) Same as Fig. 01 but for the QPO frequencies observed in SGR 
1900+14. 

Table 3: Same as Tabled but for the QPO frequencies observed in SGR 1900+14. 


QPO frequency (Hz) 

£ 

o4 e) ( Hz ) 

relative error (%) 

28 

4 

27.29 

2.55 

54 

8 

53.80 

0.36 

84 

13 

86.26 

-2.69 


Hz can be identified as £ = 4, 8, and 13 even for the calculations from the 
modified shear modulus (131) . In Fig. [6] we show such an identification for the 
stellar models with M = 1.4M 0 and R = 12 km, together with the optimal 
value of L to explain the observed QPO frequencies, i.e., L = 110 MeV. 
We also show the comparison between the observed QPO and calculated 
frequencies for the stellar models with M = 1.4M 0 and R — 12 km in 
Table [3j Additionally, we find that the observed QPO frequencies can be 
explained by the same identification even for various stellar models, where the 
corresponding optimal value of L to explain the observed QPO frequencies 
in SGR 1900+14 is shifted as shown in Fig. [0 

Although the value of L still depends on many uncertainties, it is rea¬ 
sonable to try to simultaneously explain the QPO frequencies observed in 
SGR 1806—20 and SGR 1900+14 for a specific value of L. As shown in Fig. 
El therefore, we can constrain L as 97 < L < 127 MeV by assuming that 
the central objects in SGR 1806—20 and SGR 1900+14 are neutron stars 
whose mass and radius are in the range of 1.4 < M/M 0 < 1.8 and 10 km 


< R < 14 km. Comparing to the constraint on L obtained in Ref. 31] from 


the shear modulus [Eq. (121)]. the current constraint shifts to lower values by 
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Figure 7: (Color online) Same as Fig. El but for the QPO frequencies observed in SGR 
1900+14. 


about 4 MeV due to the effect of electron screening. We also remark that 
the obtained constraint is still consistent with that from an X-ray busting 


neutron star 43 


Furthermore, we proceed to the QPO frequencies observed in SGR J1550—5418, 
which give us information independent of SGR 1806—20 and SGR 1900+14. 
Unfortunately, there are only two QPO frequencies observed, 93 and 127 Hz 
[lj. The observed frequencies, which are relatively high, might arise from 
the same oscillation mechanism other than the crustal torsional oscillations, 
but, as we shall see, the crustal torsional oscillations have little difficulty in 
reproducing the observed QPO frequencies. 

In Fig. O the frequencies of the fundamental torsional oscillations with 
l = 11 — 25 are shown as a function of L for the stellar models with 
M = 1.4M 0 and R — 12 km, together with the observed QPO frequen- 
From this figure, one finds various combinations of o tg that can re- 


cies. 


produce the observed QPO frequencies. In practice, if the observed QPO 
frequencies 93 Hz and 127 Hz are regarded as the frequencies of the funda¬ 
mental torsional oscillations with G and G(> G)> one can find that (G, G) = 
(11,15), (12,16), (13,18), (14,19), (15, 21), (16, 22), (17,23), and (18,25) are 
relevant as shown by the vertical solid lines in Fig. |9l where the correspond¬ 
ing optimal values of L are 81, 90, 101, 109, 121, 130, 139, and 161 MeV, 
respectively. Additionally, in Table U we show the possible combinations 
of and Gj the frequencies of the corresponding torsional oscillations that 
can be obtained from Eq. ([5]) by substituting the optimal values into L for 
M = 1.4M 0 and R — 12 km, and their relative errors from the observed 
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Figure 8: (Color online) The allowed values of L (shaded region) that can simultaneously 
explain the QPO frequencies observed in SGR 1806—20 and 1900+14, if the corresponding 
central objects are neutron stars having mass and radius in the range of 1.4 < M/Mq < 1.8 
and 10 km < R < 14 km. 

QPO frequencies. Considering the accuracy of the fitting formula (JE]) as well 
as the limited number of the observed QPO frequencies, we have difficulty 
in judging which combination is in best agreement with the observed QPO 
frequencies. 

Assuming that the central object in SGR J1550—5418 is a neutron star 
within 1.4 < M/M Q < 1.8 and 10 km < R < 14 km, we show in Fig. [TO] 
that several combinations of (A, 0) are consistent with the constraint on L 
from the QPO frequencies observed in SGR 1806—20 and SGR 1900+14 as 
described in Fig. [8] As in the case of SGR 1806—20 and SGR 1900+14 shown 
in Figs. [5] and [?l the upper and lower limits of the constraint on L for each 
(0,4) hi Fig. [TO] are determined by the stellar models with (M, i?) = (1.4M 0 , 
10 km) and (1.8M 0 , 14 km), respectively. From this figure, we fold that at 
least the identification of the QPOs observed in SGR J1550—5418 as the fun¬ 
damental torsional oscillations with £ = 11 and 15 is not consistent with the 
constrained region of L from the QPO frequencies observed in SGR 1806—20 
and SGR 1900+14. Additional discoveries of the QPO frequencies in SGR 
J1550—5418 would enable us to tell which of the remaining possibilities is 
best, which in turn leads to severer constraint on L. 

5. Conclusion 

In summary, we show that the observed low-lying QPOs, including the 
new ones from SGR 1806—20 and SGR J1550—5418, can be still identified as 
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Figure 9: (Color online) The frequencies of the fundamental torsional oscillations with 
£= 11 — 25, which are evaluated as a function of L for the stellar models with M = IAMq 
and R = 12 km. The horizontal broken lines denote the QPO frequencies observed in SGR 
J1550—5418. The vertical lines denote various possible correspondences between the QPO 
frequencies and the frequencies of the fundamental torsional oscillations with the values 
of £ represented by (*, *) in the figure. 


Table 4: Same as Table [2j but for the QPO frequencies observed in SGR J1550—5418, 
where various combinations of £ to explain the observed QPO frequencies can be consid¬ 
ered. _ 


QPO frequency (Hz) 

i 

o4 e) ( Hz ) 

relative error (%) 

93 

11 

93.19 

-0.21 

127 

15 

126.1 

0.73 

93 

12 

94.26 

-1.35 

127 

16 

124.8 

1.74 

93 

13 

92.66 

0.37 

127 

18 

127.4 

-0.28 

93 

14 

93.59 

-0.63 

127 

19 

126.1 

0.70 

93 

15 

91.94 

1.14 

127 

21 

127.7 

-0.57 

93 

16 

92.91 

0.098 

127 

22 

126.8 

0.15 

93 

17 

93.77 

-0.83 

127 

23 

125.9 

0.84 

93 

18 

92.61 

0.42 

127 

25 

127.5 

-0.40 
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Figure 10: The allowed regions of L to explain the QPO frequencies observed in SGR 
J1550—5418 in terms of the fundamental torsional oscillations for the stellar models with 
the mass and radius in the range of 1.4 < M/Mq < 1.8 and 10 km < R < 14 km. The 
horizontal axis l\ denotes the angular index as which the lower QPO frequency 93 Hz is 
identified, as shown in Table [I] The shaded region is the constraint on L from the QPO 
frequencies observed in SGR 1806—20 and SGR 1990+14 as described in Fig. [8] We set 
the upper limit of the vertical axis at L = 170 MeV, above which extrapolations based on 
Eq. (0 are no longer effective. 


the crustal fundamental torsional oscillations with different Us. This result 
gives a strong indication that the QPOs and the crustal modes are more or 
less related. However, there remain many problems. To make better esti¬ 
mates of the frequencies of the crustal modes, which are given as a function 
of the EOS parameter L and the neutron star mass M and radius R in the 
present work, it is indispensable to take into account effects of magnetic fields 
on the effective shear modulus. Unfortunately, it would be a tall task given 
poorly known magnetic structure. Even within the framework of purely elas¬ 
tic shear modes as considered in the present work, how to average the locally 
anisotropic shear modulus over directions for a polycrystal could modify the 
effective shear modulus significantly jlij . Also, stability of the bcc struc¬ 
ture as assumed here is endangered by fluctuations in the density of dripped 
neutrons; the resultant change in the lattice structure could affect the ef¬ 
fective shear modulus and modifications of the enthalpy density by neutron 
superfluidity j4f|. Note that with all those theoretical uncertainties, one can 
definitely assign the angular index £ to each QPO except the two QPOs ob¬ 
served in SGR J1550—5418; assignment of these two QPOs would be possible 
only after possible observations of additional low-lying QPOs from the same 
source. 
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